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Abstract. Drifting subpulse patterns in pulsar signals are frequently interpreted in terms of a model in which a 
rotating ring of sparks on the polar cap gives rise to emission from regions of the magnetsophere connected to 
the sparks by dipolar magnetic field lines. The spacing of drift-bands in time depends on the circulation rate of 
the polar cap pattern, but to first order the longitudinal phase dependence of the subpulse modulation should 
obey a frequency-independent relation determined by the geometrical configuration in a similar manner to the 
polarization position angle. We present here observations at 272-1380 MHz of PSR B0320+39 and PSR B0809-I-74, 
both of which show nearly linear drift in two longitude regions, separated by a region of reduced modula tion and 
accom panied by a large step in the phase of the subpulse pattern. We show that the observation of Bart el et alJ 
Jl98J) that the subpulse spacing for PSR B0809+74 was 1.8 times greater at 102.5 MHz than at 1720 MHz is most 
likely an artifact of the phase step, which is only present at high frequencies. The phase steps can be understood 
as a consequence of observing overlapping offset images of the polar cap spark pattern. We also detected more 
complicated, frequency-dependent behaviour that would require that the images do not simply differ by rotation 
about their centers. Detailed modelling of non-axisymmetric refraction or distorted magnetic fields is suggested 
as a means of pursuing an explanation for this phenomenon. 
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1. Introduction 

It is generally accepted that the radio emission from pul- 
sars is beamed from its points of origin along local tan- 
gents to the dipolar magnetic field lines, and that the gen- 
eral broadening of pulsar average profiles with decreas- 
ing observing frequency can be understood as a conse- 
quence of the divergence of the field lines and the scal- 
ing of emission frequency with radial distance from the 
polar cap (e.g. lKomesarofjll97(]t ICordeslll978|) . So fun- 
damental is this so-called radius-to-frequency mapping 
(RFM) to the usual ways of thinking about pulsar emis- 
sion, that many authors appear to assume that it will 
apply to longitudinal subpulse spacing (P2) in the same 
way that it applies to c omponent separations and av- 
erage profile width (e . g. iRuderm an fc Sutherland' 'l975"; 



ing subpulse patterns as manifestations of a grazing pass 
along the edge of a beam consisting of a ring ( "carousel" ) 
of beamlets circulating abou t the magnetic axis (e.g. 
iRuderman k Sutherlandlll975|) . 



Wolszczan et all Il98lt iBartel et ah. 1981; .Bartel .198h 
Izvekova et al.ll993l) . As noted bv lGil fc Krawczvkll|l996^ . 

this is only the case when the line of sight passes di- 
rectly over the magnetic pole, and the weaker depen- 
dence on frequency seen in P9 versus the ave rage pro- 
file width reported bv llzvekova et alJ l)l993(l and lGil et al.l 



ll2002t) is just as expected under the usual model of drift- 
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lEdwards fc StappersI l|2002() (ES) showed that if one 
considers the subpulse pattern as a periodic modulation 
within any given longitude interval, the phase of the mod- 
ulation as a function of longitude should be independent 
of frequency, given a circular pattern of simply shaped 
sparks (i.e. extended in either magnetic azimuth or lat- 
itude but not in both), and neglecting the effects of ro- 
tation. This is a consequence of the fact that the trans- 
formation between the polar cap pattern and the beam 
pattern as radiated from a certain height is a simple scal- 
ing of the polar opening angle (see Fig. ^ top left). The 
sampling of the polar cap effected at a particular longi- 
tude value at different radio frequencies therefore differs 
in magnetic latitude but not in magnetic azimuth, and 
for a simple circular carousel this means that the am- 
plitude of the subpulscs may differ but their phase does 
not. Figure n (bottom right) shows the difference between 
points of zero subpulse phase and points of peak subpulse 
intensity (between which P2 is measured) for a pair of 
sight-lines on a sample polar cap configuration. It shows 
that the variation with frequency of P2 is entirely a con- 
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sequence of the different amphtude windowing caused by 
the expansion or contraction of the beam. In contrast, the 
subpulse phase envelope is unaffected by the different am- 
plitude windowing and its longitude dependence should 
be describable by a simple geometric formula, indepen- 
dent of frequency (ES). The possibility therefore exists 
for a strong test of the carousel model by examining the 
longitude- and frequency-dependence of subpulse phase 
in much the same wa y as the now well-accepted geom etric 
polarization model of lRadhakrishnan fc Cook3 ( 19691) was 
tested. 

Application of phase-based analysis to PSR B0320+39 
has shown that at 328 MHz, subpulse drifting occurs in 
two distinct longitude inter vals, between which th ere is a 
phase offset of nearly 180° l|Edwards et al J 12003 ^ ESvL). 
The transition between the regions occurs too suddenly 
to be fit by the geometric model, and ESvL suggested the 
phase step reflected a shift in dominance between two su- 
perposed, phase-offset images of an underlying rotating 
subbeam system. If this is the case, in the absence of fur- 
ther complicating factors, the underlying phase envelopes 
of the two images should exhibit the predicted frequency 
independence, and any variations in the measured enve- 
lope must be due to frequency dependence of the relative 
amplitudes of the images at any given longitude. If the off- 
set is close to 180°at all longitudes, the phase must always 
be close to that of the dominant component (at that longi- 
tude) and after accounting for phase steps, the prediction 
of frequency independence can be tested. 

Similar behaviour, involving a subpulse phase jump 
accompanied by reduced modulation in the vicin- 
ity of the jump, has been r eported at 1720 and 
1420 MHz for PSR B0809-h7 4 ijWolszczan et alJ Il98lt 
IProszvnski fc WolszczanI ri986|) . although the quality of 
these results is insufRcient to strongly constrain the sharp- 
ness of the jump, or to probe for other non-linearities and 
their variation with frequency. Without this information it 
is difficult to distinguish betwee n the multiple imaging sce- 
nario and other possibilites Ce.g. lBarteligSltrDavies et al.l 
Il984ll . To investigate this phenomenon further we have 
studied PSR B0320-h39 and PSR B0809-I-74 at frequen- 
cies between 270 and 1380 MHz. 

2. Observations and Analysis 

We observed PSRs B0320-f 39 and B0809-H74 using the 
Westerbork Synthesis Radio Telescope and its pulsar back 
end, PuMa. The signals from linearly polarized receptors 
from 14 25-m dishes were decimated into bands of 10 MHz 
width, added (per-band, per-polarization) according to 
previously determined phase and amplitude factors, and 
fed in to PuMa operating in Mode 1, which acts a s a dig- 
ital filterbank; for details see IVoute et al.l l|2002[) . In all 
cases we used a configuration such that both the disper- 
sion smearing within each frequency channel and the out- 
put sampling interval were less than or equal to 409.6 /is. 
In offline processing total intensity samples from all fre- 
quency channels were combined according to the 
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Fig. 1. (Top left) Schematic diagram depicting the ori- 
gin of simultaneously observable rays (wavy lines) emit- 
ted at different altitudes. The thick lines show the path of 
plasma flows from the polar cap excitation pattern, along 
field lines to the points of emission of these rays. Due to 
REM and the divergence of the field lines, low frequency 
observers sample a more poleward path on the polar cap. 
The thin lines show other co-planar field lines, and the 
stellar surface. (Bottom right) Diagram of a carousel sys- 
tem on the polar cap. The solid circles represent a system 
of sparks circulating about the magnetic pole. The thick 
horizontal lines depict the path sampled by observers re- 
ceiving radiation emitted from two different heights in the 
magnetosphere (the upper line corresponding to a greater 
altitude) along a circular path effected by the rotation of 
the star under a specific viewing geometry. These lines are 
connected at several points to indicate simultaneously ob- 
served locations (neglecting propagation delays), spaced 
at equal time/longitude intervals At. The lines emanating 
from the center of the carousel indicate two of the twelve 
lines of zero phase in the azimuthal modulation, which are 
seen at the same pulse longitudes (±At, chosen for con- 
venience) for both observers. The small filled circles show 
the longitudes at which subpulse peaks are recorded by 
each observer (for this carousel orientation). Unlike fidu- 
cial phase points, the longitudinal separation between sub- 
pulse peaks increases with the altitude of emission. The 
overall pulse width, as indicated by the intersections of 
the outer dotted circle with the thick lines, evolves even 
more strongly. 
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appropriate compensation for interstellar dispersion, 
and where necessary adjacent samples were added to im- 
prove the signal-to-noise ratio. Based on the topocentric 
period predicted using the TEMPO^ software package and 
published ephemerides, the resultant time-series was re- 
arranged into a two-dimensional array of pulse longitude 
and pulse number for further analysis. We made obser- 
vations with bands centered at 272-302, 328, 382 and 
1345-1415 MHz, however for both pulsars, the results were 
consistent within their respective uncertainties for all fre- 
quencies between 272 and 382 MHz, so in addition to the 
1380 MHz data, we will present results only from the most 
sensitive band, which was 328 MHz in both cases. 

For pulsar signals that are dominated by a regular 
drifting subpulse pattern, the method of ES offers the best 
possible signal-to-noise ratio for the inference of the time- 
and longitude-dependence of subpulse phase. However, 
we find that for both pulsars, at 1380 MHz the regular 
drifting subpulses comprise a relatively minor portion of 
the total emission, and the time/longitude decomposition 
fails for this reason. To overcome this, in this work we 
used an alternative algorithm to ascertain the longitude- 
dependence of subpulse amplitude and phase (hereafter, 
the complex "envelope"). Specifically, we divided the 
pulse series in to 256-pulse blocks and computed Discrete 
Fourier Transforms along constant-longitude columns of 
each block to form a series of complex Longitude- Resolved 
Fluctuation Spectra (LRFS) . If the subpulse pattern is at 
least strong and coherent enough to dominate over other 
emission in its peak spectral frequency bin, that row of any 
given spectrum should differ from that of any other spec- 
trum by a simple complex factor. The phases of these fac- 
tors were determined by cross-correlation with a template 
envelope, a new template was formed by adding the peak 
rows using the measured phases, and the process repeated 
until convergence. This is essentially a coarse-grained ver- 
sion of the algorithm of ES, and has been shown to pro- 
duce consistent results in practice (ESvL) for the phase 
envelope. 

The subpulse amplitude envelope derived with this 
method is typically underestimated since the single 
Fourier coefficient is not matched to the true response 
of the signal. To overcome this, we measured the LRF 
power spectrum as the difference between on-pulse power 
spectra and off-pulse (noise) spectra, and estimated the 
variance of the subpulse-modulated part of the signal as 
the sum of power spectral values in a narrow range of fre- 
quencies around the peak response (Parseval's Theorem). 
Assuming the subpulse signal has a similar amplitude dis- 
tribution to that of a sinusoid, we then took the value 
(2i;)^/^ where v is the variance as an estimate of the 
typical height of subpulse peaks above the mean in the 
given longitude bin. As expected, for low-frequency ob- 
servations (where the subpulse signal dominates) the full 
time/longitude decomposition method (ES) gave a result 
consistent with this method. 

^ http:/ /pulsar. princeton.edu/tempo/ 



3. Results and Discussion 

3.1. PSR B0320+39 

The average profiles and subpulse amplitude and phase en- 
velopes for PSR B0320-h39 at 328 MHz and 1380 MHz are 
shown in Fig. |2 We used the same observation as ESvL at 
328 MHz, and the results are consistent within the noise 
(despite using different techniques). A sharp jump in sub- 
pulse phase is accompanied by an almost complete atten- 
uation of subpulse amplitude, consistent with destructive 
interference between two offset drift patterns with differ- 
ing longitudinal amplitude dependences. The presence of 
an additional component was the motivating factor for 
examining the subpulse properties at high frequency, with 
the prediction that it should be preceded by another jump 
to align it with the extrapolation of the subpulse phase en- 
velope in the leading part of the profile. The extra compo- 
nent is clearly detected here, as is its subpulse modulation 
(Fig[211eft, longitude ~ 26°). However, the phase envelope 
differs from the prediction, and tends towards almost sta- 
tionary modulation in the trailing component. Moreover, 
the sharp jump seen at 328 MHz (longitude ~ 19°) is ap- 
parently reduced to a smoother transition of the opposite 
sense (seen most clearly in the middle panel). The enve- 
lope resumes approximately the same phase evolution as 
the 328 MHz profile after the jump, until it is disturbed by 
the extra component (longitude ^ 23°). Clearly, the situ- 
ation is not as simple as the superposition of two rotating 
subbeam systems offset purely in magnetic azimuth, as 
earlier suggested. 

3.2. PSR B0809+74 

The average profiles and subpulse amplitude and phase 
envelopes for PSR B0809+74 at 328 MHz and 1380 MHz 
are shown in Fig. El Beginning with the 328 MHz re- 
sult, we see that the subpulse phase envelope is not as 
steep at the edges compared to other longitudes, consis- 
tent with the expectations of the carousel model. However, 
the high sensitivity of the data reveals that the phase 
behaviour is more complicated than the smooth cubic- 
like forms allowed on basic geometric grounds (ES). The 
1380 MHz phase envelope shows similar variations in sub- 
pulse phase in the leading part of the profile, but exhibits 
a striking jump in subpulse phase in the middle of the 
profile, confirming earlier resul t s of lower longitude res- 
olution (* Wolszczan et alJ Il98l[ IProszviiski fc WolszczanI 
.1986 ) . The shape of phase envelope in the trailing part 
of the profile also approximately matches that seen at low 
frequencies, with an offset of 120° in relative phase com- 
pared to the leading part. Both the large phase jump, and 
the smaller jump or steepening around 52° longitude cor- 
respond to local minima in the subpulse amplitude en- 
velope, suggesting in our view a similar interpretation to 
the phase offset in PSR B032G+39: the presence of su- 
perposed, out-of-phase drift patterns with differing longi- 
tudinal amplitude dependences. The cumulative effect of 
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both phase jumps is seen to be approximately 180°, and 
could be explained as the consequence of receiving rays 
from opposite sides of the magnetic axis, as invoked for 
PSR B0329+39 (ESvL). However, in this case the but the 
presence an additional "image" of the carousel pattern, 
offset by some amount other than 180°, is required to ex- 
plain the peak in subpulse amplitude around 55° longi- 
tude, and intermediate phase seen in the region. We con- 
firmed numerically that the major features of the complex 
subpulse envelope can reproduced using three overlapping 
Gaussian components with relative phases of 0, 35 and 
180°. However, we note that this should only be taken as 
an indication of the general suitability of this picture, and 
not a solution for this specific set of relative phase values, 
for there is considerable room for different relative phases 
given different component shapes and/or a non-linear un- 
derlying phase function under various potential viewing 
geometries. 

These results are consistent with those of previous 
studies, although for reasons outlined in the Introduction, 
w e reach somewhat d ifferent conclusions. A major result 
of iBartel et al. 1 lll98ll) was that the overall drift rate dif- 
fered by a factor ^ 1.8 between 102 and 1720 MHz, a fact 
which they used to argue that the entire profile should 
be scaled in longitude by this factor to correct for RFM. 
Performing this scaling results in a significantly wider pro- 
file for 1720 MHz than 102 MHz, on which basis it was ar- 
gued that some e mission was "missing" a t 102 MHz. For 
reasons noted by iGil fc Krawczvkl l|l996fl and reiterated 
in the Introduction to this work, RFM under the carousel 
model is not expected to alter the observed P2 value the by 
sa me factor as the pro file width, so the scaling performed 
bv lBartel etHI l)l98lh is not valid. Indeed, the value of P2 
is expected to vary relatively little, making the observed 
factor of 1.8 somewhat unexpected. However, the phase 
jump in the middle of the profile at 1720 MHz provides the 
obvious explanation for this fact. Specifically, most pairs 
of subpulses in a given pulse will lie on opposite sides 
of the phase jump and will thus appear closer together 
than expected without the phase jump. Any method of 
measurement relying on the longitudinal separation of 
subpulse peaks, in cluding the autocorrelation method of 
iBartel et al ] ill98ll) . will be subject to this distortion. The 
dotted hne in the bottoin pane l of |2 depicts the drift rate 
inferred bv lBartel etH\ l)l98l'). which can clearly be seen 
to represent a kind of interpolation between the phases at 
the peaks of the offset components we find at 1380 MHz. 
The drift rate measured at 102 MHz would correspond 
to a horizontal line in this plot, confirming that there 
is no major evolution in the drift rate between 102 and 
328 MHz. We therefore argue that the factor of 1.8 re- 
ported in subpulse spacing between 102 and 1720 MHz is 
almost entirely due to the effect of the phase offset and 
amplitude windowing. 

So me apparently anomalous results of iDavies et al.l 
l)l984|) als o receive explanation from factors discussed in 
this work. IPavies et alJ l)l984l) state that their observa- 
tion that the drift rate at 408 MHz is 1.5 times faster 



near the profile edge edge is unexpected, and may be 
caused by a non-dipole magnetic field. In fact, under the 
carousel model faster drifting is expected at the profile 
edges (ES; see also Fig. 1 this work). We suspect that the 
confusion arose due to the reciprocal relationship between 
"drift rate" (i.e. longitude displacement per time interval) 
and the longitudinal rate of change of subpulse phase (i.e. 
phase or time displacement per longitude interval). We 
alsonote that the drift speed measured by IPavies et alJ 
l)l984l) at 408 MHz actually increases monotonically with 
longitude, the asymmetry of which is not consistent with 
the model, with or without the noted misconception. In 
any case, we have shown here that the phase envelope is 
considerably more complicated than this, and yet the re- 
duced slopes at the profile edges indicate potential consis- 
tency with the carousel model if offse t multiple imaging is 
allowed. The monotonic trend seen bv lDavies et al.ll|l984^ 
probably has more to do with the perturbations between 
subpulse peaks and true points of zero phase imposed by 
the overall amplitude windowing than with the underlying 
phase variations. 

Another result of iDavies et alJ is that the sub- 

pulses at 1412 MHz apparently arrive 10-18 ms ear- 
lier than those at lower frequencies (after compensating 
for interstellar dispersion). Comparison of their Figures 
2 and 4 with Fig. [21 confirms that the offset subpulses 
they recorded appeared where our proposed alignment in- 
dicates a 120° phase offset. This correspo nds to a ver- 
tical shift of ^ 3.7 periods in Fig. 4 of IPavies et alJ 
l)l984l) . consistent with their measurement between 1412 
and 406 MHz, made with absolute time alignment. We 
may therefore conclude that our phase alignment is cor- 
rect and that high-frequency subpulses in the leading part 
of the profile are not offset from those at low frequencies. 

3.3. Implications for Models 

The subpulse phase envelopes measured for PSR 
B0320-K39 and PSR B0809-h74 are clearly not consistent 
with the smooth, symmetrical curves predicted by sim- 
ple rotating subbeam models. Moreover, they are seen 
to evolve strongly with frequency, and to exhibit sharp 
"jumps" in subpulse phase. 

Further evidence for departure from the predictions 
of a simple carousel mod e l also exists in the litera- 
ture. iHankins fc WolszczaiJ l|l987l) found support for the 
carousel model in their detection of systematic subpulse 
drift in three pulsars with "triple" profile morphologies. 
However, the subpulse phase envelopes presented for PSR 
B1918+19 clearly show a reversal in the sense of sub- 
pulse drift in the trailing component. Under the carousel 
model, points of sense reversal correspond to pulse longi- 
tudes where the line of sight is at its furthest excursions in 
the magnetic azimuth coordinate. For this to occur within 
the emission beam requires the angle between the spin axis 
and the observer to be comparable to the opening angle of 
the beam, which would result in a very broad pulse profile 
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15 20 25 30 30 40 50 60 70 80 90 

Pulse longitude (deg) Pulse longitude (deg) 

Fig. 2. Average profiles and longitude envelopes for PSR B0320+39 (left) and PSR B0809+74 (right) at 328 MHz 
(dashed lines / white circles) and 1380 MHz (solid lines / black circles). The top panel shows the average profile and 
subpulse amplitude envelope (the latter usually being lower than the average profile). The middle panel shows the 
subpulse phase envelope, plotted repeatedly with a spacing of 360°to show the inherent ambiguity. Adjacent significant 
points are joined by the shortest path. The dotted line shows the nominal phase slope (60°/° for PSR B0320+39, 
25°/° for PSR B0809-I-74). The difference between the phase envelopes and this line is plotted in the bottom panel (for 
PSR B0809-I-74, twice, with an offset of 120° to depict the rough consistency between the frequencies after correction 
for an offset ) . The dotted line in the bottom right panel indicates the slope corresponding to the P2 measured by 
iBartel et ahl l)l98j) at 1720 MHz. The error bars for PSR B0320+39 at 1380 MHz are twice the formal value. Error 
bars are omitted on other points since they are in most places very small. 
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(~ 180°). This is not the case for PSR B1918+19. Indeed, 
given the narrow width of pulsar beams, unless there is a 
trend toward alignment of spin axes and magnetic axes, 
such configurations should be very rare in the observable 
population, and sense reversals in general can be taken 
as evidence against the presence of a simple rotating sub- 
beam system. 

The clear cases of step-like features accompanied by 
reduced modulation amplitude, seen in PSR B0320-I-39 
and in PSR B0809-I-74 at high frequencies, are sugges- 
tive of the presence of superposed subbeam patterns that 
rotate together (e.g. as "images" of the same spark sys- 
tem) but are offset in the azimuthal coordinate (ESvL). 
We believe that a similar effect can explain the less or- 
dered drift behaviour seen at other frequencies in these 
pulsars, and most likely that of PSR B1918-fl9 also. For 
this to be the case, the phase envelopes of the component 
images cannot be related by a constant phase offset, re- 
quiring the corresponding subbeam patterns to have off- 
set centroids, or different (non-circular) shapes. To the 
extent that the subbeam motions deviate from concentric 
circular motion in this way, the us e of subpulse patterns 
to de termine geometric parameters llDeshpande &: RankinI 
I1999I: ES) or to "map" the po lar cap excitation pat- 
tern I Deshpande fc ri,ankinlll99fil) is not possible until the 
transformations effected in the images are determined, or 
the phase envelope is shown to be consistent with a single 
rotating pattern. 

An obvious mechanism to which to appeal is the trans- 
formation to the observer's frame, which includes aberra- 
tion and retardation. The approximate effect of these is 
to differentially shift the beam patterns in the pulse lon- 
gitude coordinate. However, as noted by ESvL, this re- 
quires that the two images be emitted at points separated 
by > 4500 km, in conflict with the current consensus con- 
cerning emission altitudes. Likewise, an altitude difference 
of > 2500 km is needed to explain the 120° phase offset 
seen in here PSR B0809-I-74 using the combined effects of 
retardation and aberration. 

Another possible explanation of the deviations ob- 

served is t hat the magnetic field is not dipolar. 

IPavies et al.l l)l984[l proposed such a condition to ex- 
plain the frequency evolution of the pulse profile of PSR 
B0809-I-74, and the offset nature of subpulses at high fre- 
quencies. Clearly this simple picture cannot now suffice, 
since a combination of offset and not offset subpulses ap- 
pears to be present. Perhaps, by attributing the multi- 
ple images to multiple discrete emission heights, some 
progress might be made with this model. 

Allowing for refraction eases the problems by opening 
the possibility of a secondary beam due to rays that cross 
the magnetic axis (ESvL). However, to produce offsets 
other than 180°, symmetry-breaking effects are needed. In 
this sense, three-dimensional refraction calculations with 
full treatment of the effects of rotation on the plasma den- 
sity and flux distribution, and retardation and aberration, 
should be seen as an essential next step in the understand- 
ing of pulsar magnetospheres. 
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